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ABSTRACT 


High  range  resolution  radar  systems  have  many  advantages  such  as  target 
classification,  resolution  of  multiple  targets,  accurate  range  measurement,  target  range 
profile  and  detection  of  low  radar  cross  section  (RCS)  targets  in  clutter.  High  range 
resolution  requires  large  bandwidths.  Stepped  frequency  waveforms  can  achieve  high 
range  resolution  by  increasing  the  effective  bandwidth  without  increasing  the 
instantaneous  bandwidth  which  would  increase  the  hardware  requirements  including 
higher  analog  to  digital  (A/D)  sampling  rates  which  are  limited  by  existing  technology. 
Under  today’s  hardware  limitations,  the  stepped  frequency  waveform  becomes  very 
important.  This  thesis  briefly  discusses  the  stepped  frequency  radar  and  associated  signal 
processing,  it  investigates  the  ambiguity  function  of  the  stepped  frequency  waveform  and 
the  stepped  frequency  radar  system.  Mathematical  expressions  of  ambiguity  functions  are 
derived  and  the  improvement  of  clutter  suppression  capability  for  the  stepped  frequency 
radar  by  rejecting  initial  pulses  is  also  discussed. 
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I.  INTRODUCTION 


High  range  resolution  radar  systems  have  many  advantages  such  as  target 
classification,  resolution  of  multiple  targets,  accurate  range  measurement,  target  range 
profile  and  detection  of  low  radar  cross  section  (RCS)  targets  in  clutter.  High  range 
resolution  is  typically  achieved  either  by  decreasing  the  duration  of  the  transmitted  signal 
or  by  modulating  a  transmitted  signal  of  relatively  longer  duration.  In  either  case,  the 
instantaneous  bandwidth  of  the  signal  goes  up  and  thereby  increases  the  hardware 
requirements  to  include  higher  A/D  sampling  rates  (which  are  limited  by  existing 
technology).  However,  with  the  stepped  frequency  waveform,  it  is  possible  to  achieve 
higher  range  resolutions  with  lower  instantaneous  bandwidths  and  lower  A/D  sampling 
rates.  Irrespective  of  the  waveform  or  compression  method  used,  high  range  resolution 
requires  large  bandwidths.  For  the  stepped  frequency  waveform,  a  large  bandwidth  is 
obtained  sequentially  by  changing  the  carrier  frequency  in  steps  over  several  pulses  instead 
of  within  a  single  pulse.  This  waveform  potentially  can  be  implemented  in  high 
performance  radar  systems  with  lower  A/D  sampling  rates,  which  is  one  of  the  current 
bottlenecks  in  the  development  of  radar  systems.  Because  of  the  importance  of  the 
stepped  frequency  waveform,  this  thesis  briefly  discusses  the  stepped  frequency  radar,  its 
associated  signal  processing,  and  analyzes  it  by  the  ambiguity  function.  In  Chapter  II  the 
stepped  frequency  waveform  and  its  corresponding  radar  system  are  introduced.  In 
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Chapter  in  the  matched  filter  and  the  correlation  definition  of  the  ambiguity  function  are 
first  discussed.  The  properties  and  applications  of  the  ambiguity  function  are  also  briefly 
discussed.  The  mathematical  expression  of  the  auto-ambiguity  function  for  the  stepped 
frequency  waveform  is  then  derived  from  the  definition  and  is  verified  by  comparing  the 
ambiguity  diagram  with  one  obtained  by  simulation.  Finally,  the  mathematical  expression 
of  the  cross-ambiguity  function  for  the  stepped  frequency  radar  signal  processor  is  derived 
and  verified,  and  the  improvement  of  the  clutter  suppression  capability  for  the  stepped 
frequency  radar  by  the  rejection  of  initial  pulses  is  discussed.  In  Appendix  A,  the 
ambiguity  functions  for  the  single  pulse,  constant  frequency  pulse  train,  linear  frequency 
modulated  (LFM)  pulse,  and  the  discrete  frequency  modulated  pulse  are  given  for  the 
purpose  of  comparison.  Mathematical  expressions  of  the  ambiguity  functions  for  these 
waveforms  are  derived  and  verified.  In  Appendix  B,  the  computer  programs  of  ambiguity 
functions  for  all  the  waveforms  discussed  in  the  thesis  are  listed. 
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H.  OPERATING  PRINCIPLES  OF  THE  STEPPED  FREQUENCY  RADAR 

A.  STEPPED  FREQUENCY  WAVEFORM 

The  stepped  frequency  waveform  can  be  described  by  a  series  of  N  coherent 
pulses,  where  the  carrier  frequency  increases  from  pulse  to  pulse  by  a  fixed  increment  A f 
as  shown  in  Figure  1.  Each  pulse  has  a  fixed  pulse  width  Ts,  and  the  pulses  are  transmitted 
at  a  fixed  pulse  repetition  interval  (PRI).  The  frequency  of  the  kth  pulse  is  given  by 

fk=T0+(k-l)A/,  (1) 

where  f0  is  the  nominal  carrier  frequency  and  A/ is  the  frequency  step  size.  [Ref.  1,  pp. 
160-161] 

The  instantaneous  bandwidth  of  this  waveform  is  approximately  equal  to  the 
inverse  of  the  pulse  width  and  is  much  less  than  the  effective  bandwidth.  The  waveform's 
effective  bandwidth,  denoted  as  Beff,  is  determined  by  the  product  of  the  number  of  pulses 
N,  and  the  frequency  step  size  A/ as  given  below. 

Brff=NAf  (2) 

The  range  resolution  for  any  waveform  is  dependent  on  the  effective  bandwidth  of  the 
waveform.  The  range  resolution  of  the  stepped  frequency  waveform  is  given  as 

Ar=2^  (3) 

where  c  is  the  speed  of  the  light.  Thus  the  range  resolution  can  be  made  finer  by  increasing 

the  effective  bandwidth  NAf.  [Ref.  2,  p.  234] 
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1  2  3  4  N 

fO  fO+Af  fO+2Af  fD+3Af  fO+(N-l)Af 


N  :  number  of  pulses  Tr  :  pulse  repetition  interval  (PRI) 

Ts  :  pulse  width  Af ;  frequency  step  size 

fO  :  nominal  carrier  frequency 

Figure  1 .  The  Stepped  Frequency  Waveform 
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The  effective  bandwidth  is  achieved  by  the  signal  processing  of  the  data  obtained 
during  the  coherent  processing  interval  (CPI),  which  usually  is  the  time  on  target.  The  CPI 

is  given  by  the  product  of  the  number  of  pulses  N  and  the  pulse  repetition  interval  (PRI). 

CPI=NxPRI  (4) 

B.  SYSTEM  DESCRIPTION 

The  block  diagram  of  a  radar  which  employs  a  stepped  frequency  waveform  is 
shown  in  Figure  2.  The  transmitted  signal  in  this  radar  is  generated  as  in  conventional 
coherent  radars  by  mixing  RF  frequency  (STALO)  and  IF  frequency  (COHO).  However, 
the  frequency  step  in  each  pulse  is  added  by  the  frequency  synthesizer  as  shown  in  Figure 
2.  On  the  receiver  side,  a  reverse  operation  of  taking  out  the  RF,  step,  and  IF  frequencies 
is  performed.  The  IF  amplifier  serves  as  a  matched  filter  to  the  envelope  of  each  individual 
pulse.  The  output  of  the  synchronous  detector  is  a  video  signal  in  the  form  of  in-phase  and 
quadrature  (I-Q)  components.  The  output  of  the  A/D  is  stored  and  organized  according  to 
range  bins.  Received  signals  from  all  N  pulses  in  each  range  bin  is  collected.  This  data  may 
be  weighted  before  taking  the  DFT.  The  output  of  the  DFT  represents  a  finely  resolved 
range  bin.  The  size  of  the  original  range  bin  is  cx/2.  However,  after  DFT,  the  width  of  the 
subdivision  of  a  finely  resolved  range  bin  is  cx/2NAfx,  which  is  equivalent  to  having  a 
bandwidth  of  NAf  [Ref.2,  p.  236].  The  unambiguous  range  R^  for  this  waveform  is  c/2Af. 
This  will  be  quite  short  but  it  may  not  be  a  problem  as  it  is  applicable  within  the  same 
range  bin.  For  avoiding  wrap  around  the  target,  Ru  should  be  greater  than  the  target 
extent. 
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Figure  2.  System  block  diagram  for  the  Stepped  Frequency  Radar 
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m.  ANALYSIS  OF  THE  STEPPED  FREQUENCY  RADAR 
BY  THE  AMBIGUITY  FUNCTION 

A.  INTRODUCTION  TO  AMBIGUITY  FUNCTION 

1.  Matched  Filter 

The  matched  filter  maximizes  the  output  peak  signal  to  noise  power  ratio  for 
the  detection  of  a  known  signal.  In  the  presence  of  additive  white  Gaussian  noise,  the 
impulse  response  of  the  matched  filter  is  the  complex  conjugate  of  the  image  of  the 
received  waveform,  h(t)=s*(t1-t),  where  s(t)  is  the  received  signal,  that  is,  it  is  the  same  as 
the  received  signal  run  backward  in  time  starting  from  fixed  time  tr  [Ref.  3,  p.  261]  In 
order  to  simplify  the  mathematical  developement,  it  is  customary  to  set  tj  at  zero.  The 
output  of  the  matched  filter  is  then  written  as 

OO  OO 

y(t)  =  f  s(X)h(t- X)dk  =  J  s(X)s*  (k  -  t)dk  (5) 

— -oo  — oo 

2.  Definition  of  the  Ambiguity  Function 

If  the  transmitted  signal  St(t)  is  u(t)ej2l?ot ,  then  received  signal  Sr(t)  is  given 

by  u(t-TD ,  where  To  is  the  unknown  round  trip  time  delay  given  by  ^r- 
and  fd  is  the  Doppler  shift  of  the  returned  signal.  Assuming  that  the  radar  receiver  is 
matched  to  the  transmitted  signal,  we  obtain  the  output  of  the  matched  filter  as 

So[td)  =2  u(t  -  To  xu*(t-TD  yM‘-T>)dt  ,  (6) 

/ 

where  TD  is  the  estimate  of  time  delay. 
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/  / 

It  is  customary  to  set  T D  and  fo  equal  to  zero,  and  to  define  Td~Td  as  x.  The 

output  of  the  matched  filter  is  then 

So(t)  =  e~j2^dZx  J  u(t-x)u*(t)ej2^dtdt.  (7) 

— oo 

Since  the  interest  is  only  in  the  magnitude  of  the  matched  filter  output, 
neglecting  the  first  exponential  term,  we  obtain  the  output  of  the  matched  filter  as 

X(x ,fd)  =J  u(t- z)u* (t)e^dtdt  (8) 

— OO 

In  this  form  a  positive  x  indicates  a  target  beyond  the  reference  delay  t0 , 
and  a  positive  /d  indicates  an  incoming  target  [Ref.  4].  In  the  literature  the  term 
"Ambiguity  Function"  is  interchangeably  used  for  X(x/d),  IX(x/d)l,  and  IX(x/d)l2.  In  this 
thesis,  the  plot  of  IX(x/d)l2  is  called  the  ambiguity  diagram.  If  the  receiver  filter  is  matched 
to  the  transmitted  signal,  IX(x/d)l2  is  termed  the  auto-ambiguity  function,  otherwise,  it  is 
termed  the  cross-ambiguity  function. 

3.  Properties  of  the  Ambiguity  Function 

The  function  IX(x,/d)l2  has  the  following  properties  [Ref.  5,  p.  412]: 


Maximum  value  of  \X(x,fd)\2  =  lx(0,0)l2  =  (2 E)2, 

(9) 

\X(-x,-fd)\2  =  \X(x,fd)\2, 

(10) 

IX(t,0)I2  =  I  J  u{t-x)u*{t)dt\2 , 

— OO 

(11) 

\X(0Jd)\2  =  1  J  u2{t)e]1^ddt\2 , 

— oo 

(12) 
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J  J  \X(xJd)\2dxdfd  =  (2 E)2 .  (13) 

The  first  equation  given  above.  Equation  9,  states  that  the  maximum  value  of 

the  ambiguity  function  occurs  at  the  origin  and  its  value  is  (2E)2,  where  E  is  the  energy 
contained  in  the  echo  signal.  Equation  10  shows  that  the  ambiguity  function  is  symmetric 
about  the  origin.  Equation  11  describes  that  the  ambiguity  function  along  the  time  delay 
axis  is  the  auto-correlation  function  of  the  complex  envelope  of  the  transmitted  signal. 
Equation  12  shows  that  along  the  frequency  shift  axis  the  ambiguity  function  is 
proportional  to  the  spectrum  of  u2(t).  Equation  13  states  that  the  total  volume  under 
the  ambiguity  function  is  a  constant  equal  to  (2E)2. 

The  ideal  ambiguity  diagram  consists  of  a  single  peak  of  infinitesimal  thickness 
at  the  origin  and  is  zero  everywhere  else,  as  shown  in  Figure  3.  The  single  spike  eliminates 
any  ambiguities,  and  its  infinitesimal  thickness  at  the  origin  permits  the  frequency  and  the 
echo  delay  time  to  be  determined  simultaneously  to  as  high  a  degree  of  accuracy  as 
desired.  It  also  permits  the  resolution  of  two  targets  no  matter  how  close  together  they  are 
on  the  ambiguity  diagram.  Because  of  two  restrictions  the  ideal  ambiguity  function  is  not 
possible;  first,  the  maximum  height  of  the  ambiguity  function  must  be  equal  to  (2E)2, 
which  is  limited,  and  second,  the  total  area  under  the  ambiguity  function  must  be  finite 
and  equal  to  (2E)2.  A  reasonable  approximation  to  the  ideal  ambiguity  function  might 
appear  as  shown  in  Figure  4.  This  waveform  does  not  result  in  ambiguity  since  there  is 
only  one  peak,  but  the  single  peak  might  be  too  broad  to  satisfy  the  requirements  of 
accuracy  and  resolution.  If  the  peak  is  made  too  narrow,  the  requirement  for  a  constant 
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Figure  3.  Ideal  ambiguity  diagram  "after  Ref.  [5], 
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volume  might  cause  peaks  to  form  at  regions  of  the  ambiguity  diagram  other  than  the 
origin  and  give  rise  to  ambiguities.  Thus  the  requirements  for  accuracy  and  ambiguity  may 
not  always  be  simultaneously  possible  to  satisfy.  [Ref.  5,  pp.  412-413] 

The  particular  waveform  transmitted  by  a  radar  is  chosen  to  satisfy  the 
requirements  for  (1)  detection,  (2)  measurement  accuracy,  (3)  resolution,  (4)  ambiguity, 
and  (5)  clutter  rejection.  The  ambiguity  diagram  may  be  used  to  assess  qualitatively  how 
well  a  waveform  can  achieve  these  requirements.  The  maximum  value  of  the  ambiguity 
function  is  an  indication  of  the  radar’s  detection  capabilities.  Since  the  ambiguity  diagram 
is  often  normalized  so  that  the  maximum  value  is  equal  to  one,  the  ambiguity  diagram  is 
seldom  used  to  assess  the  detection  capabilities  of  the  waveform.  The  accuracy  with  which 
the  range  and  the  velocity  can  be  measured  by  a  particular  waveform  depends  on  the 
width  of  the  spike,  centered  at  IX(0,0)l2,  along  the  time  and  frequency  axes.  The 
resolution  is  also  related  to  the  width  of  the  central  spike,  but  in  order  to  resolve  two 
closely  located  targets  the  central  spike  must  be  isolated.  It  cannot  have  any  high  peaks 
nearby  that  might  mask  another  target  close  to  the  desired  target.  A  waveform  that  yields 
good  resolution  will  also  yield  good  accuracy,  but  the  reverse  is  not  always  true.  The 
presence  of  additional  spikes  can  lead  to  ambiguity  in  the  measurement  of  target 
parameters.  An  ambiguity  measurement  is  one  in  which  there  is  more  than  one  choice 
available  for  the  correct  value  of  a  parameter,  but  only  one  choice  is  appropriate.  Thus  the 
correct  value  is  uncertain.  The  ambiguity  diagram  may  be  used  to  determine  the  ability  of  a 
waveform  to  reject  clutter.  If  the  transmitted  waveform  is  to  have  good  clutter-rejection 
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properties,  the  ambiguity  diagram  should  have  little  or  no  response  in  the  regions  of  clutter 
[Ref.  5,  pp.  418-420].  Some  examples  of  the  ambiguity  function  are  shown  in  Appendix  A 
which  includes  single  pulse,  pulse  train,  linear  frequency  modulation  pulse,  and  stepped 
frequency  modulation  pulse. 


B.  DERIVATION  OF  THE  AUTO-AMBIGUITY  FUNCTION  FOR  THE 

STEPPED  FREQUENCY  RADAR 


In  order  to  derive  the  auto-ambiguity  function,  according  to  the  definition, 

Equation  8,  we  need  the  expression  of  the  complex  envelope  of  the  transmitted  signal.  The 

transmitted  signal  of  the  stepped  frequency  radar  consisting  of  a  train  of  N  pulses  and 

changing  carrier  frequency  from  pulse  to  pulse  by  fixed  increment  A/  is  shown  in  Figure 
1 .  It  is  mathematically  represented  as 


St(t)  =  at  X1  u(t  -  nTr)e^+nm  ?  (14) 

n=0 

where  u(t~nTr)  =  l,  if  nTr  <t<  nTr  +  Ts , 

0,  elsewhere,  (15) 

fstep  =  nAf  ,  n  =  0,  1  ,and 

fo  = fcoho  ^  f state , 

fcoho  :  the  frequency  of  the  stable  coherent  oscillator , 
f staio  :  the  frequency  of  the  stable  local  oscillator , 
fstep  :  step  frequency, 

A/  :  the  frequency  step  size, 

N  :  the  number  of  the  pulse  sequences, 

Tr  :  pulse  repetition  interval  (PRI), 

Ts  :  pulse  width, 

at  :  the  amplitude  of  the  transmitted  signal. 
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The  transmitted  signal  can  be  rewritten  as 


St(t)  —  at 


X1  u(t  -  nTr)eJ'lm^t 

.  n=0 


*Wo* 


=  atU{i)e^t  (16) 

AM 

where  U(t)  =  X  u(t  -  nTr)eJl7lnAfl  is  the  complex  envelope  of  the  pulse  sequence.  The 

n- 0 

signal  received  due  to  a  discrete  scatterer  can  be  written  as 

Sr(t)  =  arUit-xje&totfMri  (17) 

where 


ar  :  the  amplitude  constant, 

x  :  the  round  trip  time  between  the  radar  and  the  scatterer, 

fd  :  the  Doppler  frequency  shift  due  to  the  normal  relative  velocity  between  the 

radar  and  the  scatterer, 

x  =  f,  (18) 

and 

U  =  2jrh-  (19) 

Substituting  the  complex  envelope,  U(t)  ,  of  the  transmitted  signal  into 

Equation  8,  we  obtain  the  auto-ambiguity  function  of  the  stepped  frequency  radar  as 


X(x ,fd)  =  J  X1  u(t-mTr - t)^k)xY  ut(t-riTr)e~]'lml^  x e^dt.  (20) 

— oom=0  n=0 


/ 

Changing  variable  t-nTr  by  t  ,  we  change  the  above  equation  to 


N-l/V— i 

X(xjd)  =  X  X  e~j2nn2Afr r  x  ej2nmn^T  r  X  e~j2rLmAft  x  ej2^dnTr 

m=0/t=0 


x  J  u*  (t  Iwf/  -irn-n)  -  X dt  . 


(21) 
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Defining  the  round  trip  time  T  in  terms  of  integer  and  fraction  parts 


x  =  (p  +  y)Tr,  (22) 

where  p  is  an  integer,  y  is  the  fraction,  and  0  <  y  <  1 ,  we  obtain  the  ambiguity  function  as 

X(x,fd)  =  (Ts  -  iTr)  X  e^-P^Ts^Tr)  x  ej2npAfi  x 

i WrPWriTr) 


min  (N—  1 JV-1  +p) 

:  2  eP-ltotfd-p&fiTr  x  g-jlltnAft 

n= max  (0j>) 


for  0<y<£, 


(23-1) 


:  (J  +  (y-  DY  ■)  x  «W<r(p+ DAfflrKr-Wr)  X 
v  '  '  jc(^-(p+i)y)(rs+(r-i)rr) 


min(iV-l  jV+z?) 
/i=max  (Oj?+l) 


for  1  <y<  1 ,  and 


'0, 


Tr 


elsewhere. 


(23-2) 


(23-3) 


After  simplifying  the  summation  term  and  taking  the  absolute  value,  we  obtain  the 
expression  for  the  auto-ambiguity  function  as 

I V/V  r  \  I  _  /-t  „/r  \  w  I  sinWd-pWHTs-iTr)  |  w  |  sin K<fd-(2pYf>y)(N-\p\ )Tr  \ 

\X(Xjd)\  -  (Ts  -  YTr)  X I  X  I . '  ’ 


for  0<y<^f, 


(24-1) 
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=  0 ,  elsewhere .  (26-3) 

One  can  easily  show  that  the  auto-ambiguity  function  of  the  stepped  frequency 

radar  will  become  the  ambiguity  function  of  the  pulse  radar  (shown  in  Appendix  A)  by 

setting  A/  to  zero.  The  ambiguity  diagram,  that  is  the  plot  of  IX(t/d)l2  for  the  stepped 

frequency  waveform  for  the  following  parameters:  number  of  pulses  N  equals  to  4, 
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duty-cycle  equals  0.2,  and  frequency  step  size  Af  equals  2PRF  is  plotted  in  Figure  5.  This 
plot  somewhat  looks  like  the  ambiguity  function  of  the  pulse  train  rotated  at  an  angle.  The 
contour  plot  of  the  ambiguity  function  is  shown  in  Figure  6.  Plots  of  a  cut  along  the 
frequency  axis  for  t=0  and  a  cut  along  the  time  delay  axis  for/d=0  are  shown  in  Figures  7 
and  8.  It  can  be  seen  from  Figure  7  that  the  spikes  appear  at  multiples  of  PRF  in  the 
frequency  axis.  The  first  null  for  each  spike  is  l/(NTr)  away  from  the  center  of  the 
spike.  Therefore,  each  spike  has  null  to  null  width  of  2/(NTr)  along  the  frequency  axis. 
Amplitudes  of  those  spikes  vary  as  a  SINC  function.  Zero  amplitude  will  appear  at  l/Ts. 
From  Figure  8,  there  is  2N-1  spikes  along  the  time  delay  axis.  From  Eq.  26-1,  it  is  easy 
to  show  that  the  first  null  for  the  central  spike  (p=0)  occurrs  at  l/(NAf).  Accordingly,  if 
1  /(NAf)  is  less  than  pulse  width  Ts ,  the  stepped  frequency  radar  will  have  better  range 
resolution  than  the  pulse  radar.  Plots  of  cut  along  the  delay  axis  with  different  frequency 
step  Af, s  are  presented  in  Figure  9.  The  larger  the  radar's  frequency  step,  the  better  the 


resolution  of  the  radar. 


o  o 


frequency  (PRF)  '10  '4  delay  (PRI) 


Figure  5.  Auto-ambiguity  Diagram  of  the  Stepped  Frequency  Radar, 
Pulse  Number  N=4,  Duty-cycle=0.2,  A/=2PRF. 
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frequency  (PRF) 


delay  (PRI) 


Figure  6.  Contour  Plot  of  the  Auto-ambiguity  Diagram  for  the  Stepped 

Frequency  Radar,  Pulse  Number  N=4,  Duty-cycle=0.2  ,A/=2PRF. 
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-  5  0  5 

frequency  (prf) 

Figure  7.  Plot  of  Cut  along  the  Frequency  Axis  for  x  =  0  of  the  Auto-ambiguity  Function 
for  the  Stepped  Frequency  Radar,  Pulse  Number  N=4,  Duty-cycle=0.2  ,  Af= 


Figure  8.  Plot  of  Cut  along  the  Time  Axis  for  f=0  of  the  Auto-ambiguity  Function  for  the 
Stepped  Frequency  Radar,  Pulse  Number  N=4,  Duty-cycle=0.2,  A/=2PRF. 
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C.  DERIVATION  OF  THE  CROSS-AMBIGUITY  FUNCTION  FOR  THE 


STEPPED  FREQUENCY  RADAR 

In  this  section  the  mathematical  expression  of  the  ambiguity  function  for  the  stepped 
frequency  radar  receiver  and  processor  as  shown  in  Figure  2  are  derived.  Since  the  filter 
used  in  this  system  is  only  matched  to  the  individual  pulse  instead  of  matched  to  the 
transmitted  signal  which  is  composed  of  N  pulses,  the  ambiguity  function  of  this  stepped 
frequency  radar  processor  is  called  the  cross-ambiguity  function.  The  stepped  frequency 
radar  can  be  mathematically  modeled  as  shown  ih  Figure  10.  The  transmitted  and  received 
signals  were  defined  in  Equation  16  and  Equation  17  respectively.  The  received  signal  is 
mixed  with  a  local  oscillator,  stepped  frequency  synthesizer,  and  a  coherent  oscillator.  The 
stepped  frequency  synthesizer  is  synchronized  so  that  the  transmitter  and  receiver  are  on 
the  same  frequency  step.  As  a  result,  multiple  time  around  echoes  will  have  frequencies 
which  vary  by  multiples  of  frequency  step  size  A/.  The  signal  is  then  passed  through  a  filter 
which  is  matched  to  a  single  pulse  envelope.  The  output  of  the  pulse-envelope  matched 
filter  is  sampled  and  the  samples  undergo  the  discrete-time  signal  processing  shown  in 
Figure  10. 

The  output  of  the  stepped  frequency  synthesizer  used  in  the  receiver  is  defined  as 

Ss'epit)  |i(i - nTr)e-J2^1 ,  (27) 

«=0 

where 

\iit-nTr)  =  1,  if  nTr  <t<(n+ 1)77 ,  and 

0,  elsewhere.  (28) 
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The  output  of  the  matched  filter,  z(t),  can  be  derived  by  taking  the  convolution  of 
the  impulse  response  of  the  matched  filter,  h(t),  and  the  received  signal  at  the  input  of  the 
matched  filter.  The  received  signal  at  the  input  of  the  matched  filter  can  be  defined  as 

Si(t)  =  e~j2^ot  x  Sstep ( t )  x  Sr(t) .  (29) 


The  above  equation  can  be  rewritten  as 


z(to)  =ar  j  e  j2v?°x x ej2l^d(- ^  xJxh*(to~ r|)<ir| , 


where 


J  =  L  (X(ti - mrr)e-^xX  h(t|-t- nTr)eJ2m^^ 

m=0  rv= 0 


(31) 

(32) 


and  x  is  defined  in  Equation  22. 


If  each  pulse  of  the  received  signal  is  assumed  to  fall  within  the  same  pulse  of 


the  stepped  frequency  signal  no  matter  how  x  varies,  then  J  can  be  expressed  as 


min  (N-l-pfi~l) 

J  =  e-j2KKm  x  u(r\-x-nTr)e-j2Kn^t, 

n~ max  (0 -p) 


where  k  =p,  if  pTr  <x<  (p+  l)Tr-Ts, 

=p  + 1,  if  (p  +  l)7Y-r,<T<(p+l)7Y. 


Substituting  Equation  33  into  Equation  31,  we  obtain  the  output  of  the  matched 


filter  as 


z(t0)  =  are-j2^+f^x  J 


min  (A^—  1  — /?  JV—  1 ) 

X  u(T)  -x- nTr)e-j2nn^  x  A* (t0  -  T])dn . 

m=max  (0 -/?) 


Letting  E,  =  q  -  to ,  we  rewrite  the  output  of  the  matched  filter  as 


z(to)  =  are~j2n(fo+^'z  xej2'n^d~K^,t° 


00  P  mm(/V~l~/7  jV-l) 

xj  ej2K(f*-^  X  ui^  +  to-i-nTr)  xh*(-Qd^ 
-o®  L  «=max  (0 -/?)  J 


=  are  ftWo+fa)*  x  ej2^mt°x  J  x  F(£,  -(x-to),x)  x  A*  , 

—CO 

where  =/<*  -  kA/ , 


min  (AM— pjV— 1) 

K(t,x)=  X  u(t-nTr)e~J2m & . 

/i=max(0 /?) 
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It  is  customary  to  set  t0  to  0  for  the  computation  of  the  ambiguity  function  and 
neglect  the  first  exponential  term,  since  we  are  only  interested  in  the  magnitude  of  the 
matched  filter  outputs.  Redefining  the  output  of  the  matched  filter  z(t0)  as  X(x,fm),  we 

obtain  the  output  of  the  matched  filter  as 

X(z,fm )  =  J  x  V(t  -  x,  x)  x  h* ( -t)dt .  (39) 

-~oo 

According  to  the  definition  of  the  ambiguity  function  in  Equation  8,  we  treat  fm  as 
the  frequency  shift  of  signal  at  the  input  of  the  matched  filter.  This  frequency  shift  consists 
of  the  Doppler  shift  and  the  frequency  changes  due  to  the  stepped  frequency  waveform. 
V(t,x)  is  then  the  envelope  of  signal  at  the  input  of  the  matched  filter.  Thus  X(x/m) 
becomes  the  cross-ambiguity  function  of  the  stepped  frequency  radar  without  the 
processor.  In  this  form  a  positive  x  indicates  a  target  beyond  the  reference  delay  t0. 

In  the  rest  of  the  section,  the  ambiguity  function  is  extended  to  include  the  signal 
processor.  The  matched  filter  output  is  sampled  at  instants  t=n Ts  and  the  samples 
undergo  discrete-time  signal  processing.  The  sample  rate,  1/TS,  is  greater  than  or 
equal  to  the  bandwidth  of  the  pulse  envelope  u(t).  Further,  the  pulse  repetition  interval 
(PRI),  Tr ,  is  an  integer  multiple,  k,  of  Ts.  The  sample  stream  for  all  the  pulses  is  sorted 
into  k  parallel  processing  streams  termed  ambiguous  range  bins.  Processing  within  each 
range  bin  utilizes  only  samples  from  time  instants  separated  by  integer  multiples  of  Tr. 
[Ref.  6] 
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The  ith  sample  for  range  bin  k  can  be  defined  as 

zk(i)  =  z(kTs  +  iTr).  (40) 

These  samples  from  the  matched  filter  output  which  are  associated  with  a  range 

bin  are  processed  together  by  a  discrete  Fourier  transformation  (DFT)  processor.  A 

Hamming  window  is  applied  before  executing  the  DFT  to  reduce  frequency  sidelobes. 

Out  of  N  samples  in  a  range  bin,  the  first  I  samples  are  rejected  to  allow  for  the  settling  of 

transients  due  to  range  ambiguous  clutter  [Ref.  7].  DFT  is  taken  as  the  last  M  weighted 

samples  of  each  range  bin  as  given  below 

Z(k,  l)  =  L  w(i)e~j2niM zk(i +/) ,  (41) 

i= 0 

where 


k  :  range  bin  index, 
l  :  Doppler  filter  index, 

M  :  the  number  of  samples  that  DFT  processed,  and  M  <N- \p\  , 

i  :  sample  index,  i=0,l,...,M-l, 
w  :  the  window  coefficients, 

I  :  the  number  of  samples  that  will  be  rejected  initially. 


Substituting  Equations  35  and  40  into  Equation  41,  we  obtain  the  /th  output  of  the 

DFT  for  range  bin  k  as 

Z(k,l)  =  are-J2K^xx 


00  min  (AM  -p  /VM ) 

J  eW**\  2  u(x\  —  t  —  x  h*  (kTs  -r\)dT], 

—=o  n- max  (0~p) 

where 


M-l  n  .  / 

K  it)  =  I  ( t+  (i  +r)Tr) . 


(42) 


(43) 
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is  the  impulse  response  of  the  radar  receiver  processor  which  includes  the 
pulse  envelope  matched  filter,  weighting  function  and  Doppler  processor . 

Equation  42  has  the  same  form  that  Equation  35  does.  So  after  changing  variables, 
setting  kTs  equal  to  zero  and  neglecting  the  first  exponential  term,  the  /th  output  of  DFT 
can  be  written  as 

=2  e^‘V(t-x,x)h*(-t)dt,  (44) 

The  above  equation  defines  the  cross-ambiguity  function  associated  with  the 
mathematical  model  of  the  stepped  frequency  radar  processor  as  illustrated  in  Figure  10. 
The  subscript  l  indicates  that  a  distinct  ambiguity  function  is  defined  for  each  Doppler 
filter. 

For  the  received  pulse  to  be  matched  to  the  transmitted  pulse,  h(t)  should  be  equal 
to  j-i(-t).  From  Equations  37,  42,  and  43,  after  simplifying,  the  ambiguity  function  for  /th 
DFT  can  be  expressed  as 

min  (N-l-l+prN-l-IrM~l)  .  / 

Xi(x,fd)  =  Z  w(i)  x  e~j2iu^  x  e-JW+HWi  x 

i=max  (0^7-/) 


x(Ts  -  iTr)  x  e^-KA/KW,)  x 


sfcntffj-TSyyJs-iTr) 

K(fd-Kty)(Ts^/Tr) 


for  0  <  y  <  Yr , 


(45-1) 
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min  (N-I+p  ft~I  Jtf- 1 )  . 

Xiix,fd)  =  2  wii)  x  X  x 

j=max  (0,/H-l-/) 


xrrt  +  (y—  l'lTV')  x  ertfd-W)(Ts+(y-i)Tr)  x  )rr) 

v  5  u  '  r'  7t(^-K4rt(rJ+(r-i)rr)  » 


for  1  -  j-  <  y  <  1 , 


(45-2) 


elsewhere. 


(45-3) 


Taking  the  absolute  value  of  the  above  equation,  we  obtain  the  ambiguity  function 


\X,(XM\  =  (TV  -y7V)  x  I  lx 


min  (iV-l-/+p  jV- 1  -/  JVf- 1 )  _  .  * 

2  vv(z")  x  x  ejlK{fd~i2p^)Af)iTr  I , 

r=max  (0,/?-/) 

for  0  <  y  <  ^ , 


(46-1) 


(T  t-(v  nr  1  s™TC(^(p+b4m7>CH)7V)  i 

(J,  +  (Y  i)ir)XI  IX 


min  (N-I+p  ft-IJtf-l)  . , 

X  w(i)  x  x  ePWa-Qp+f+imKr  |  ^ 

z=max  (0^7+1-/) 

for  1  -j-  <y<  1 , 


(46-2) 


elsewhere. 


(46-3) 


The  above  equation  becomes  the  cross-ambiguity  function  of  the  stepped 


frequency  radar  processor.  The  cross- ambiguity  function  can  describe  the  auto-ambiguity 
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function  of  the  stepped  frequency  radar  by  setting  the  Doppler  filter  index  /  to  zero  and 
assuming  a  uniform  window,  w(i)=l,  setting  /  to  zero,  and  setting  the  number  of 
transmitted  pulses,  N,  equal  to  the  number  of  samples  processed,  M.  Under  these 
assumptions  and  after  simplifying,  the  ambiguity  function  becomes 


\Xo(x,fd)\=(Ts-iTr)x\ 


StoWd-pAflCTs-iTr) 

W<r-pAf)(Ts-*(Tr) 


ixl 


sin7C(/^2/H7)A/)(jV-l/>l)rr  j 
sin  r 


forO<y<|f, 


(47-1) 


=  (T,  +  (y-1)7V)x 


j  smTtifj-lp+rmTMy-DTr)  I  j  sin 7i(/j-(2p+7+l)A/)(A'-lj7+ll)rr  i 
Wd-(P+l)mTMy-l)Tr)  X  sin3t(/'d-(2ip+Y+l)A/)rr  1 


for  1  -  <  y  <  1 ,  (47-2) 


=  0, 


elsewhere. 


(47-3) 


The  above  equation  is  exactly  the  same  as  the  expression  of  the  auto-ambiguity 
function  in  Equation  24.  If  we  let  the  frequency  step,  A/,  be  zero,  it  is  easy  to  prove  that 
Equation  38  becomes  the  ambiguity  function  of  a  pulse  train  as  expressed  in  Appendix  A 
Equation.  A.  13. 

The  ambiguity  diagram  of  the  0th,  1=0,  Doppler  filter  of  the  stepped  frequency 
radar  with  a  Hamming  window  in  the  Doppler  processor  is  illustrated  in  Figure  11.  Figures 
12  and  13  are  plots  of  cut  along  the  frequency  axis  and  delay  axis  respectively.  From  the 
comparison  of  Figures  12  and  13  with  Figures  7  and  8  which  are  the  cut  plots  of  the 
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frequency  (prf) 


-10  -4 


delay  (pri) 


Figure  1 1.  Cross-ambiguity  Diagram  of  the  Stepped  Frequency  Radar,  Pulse  Number 

N=4,  Number  of  Pulses  Processed  M=4,  Duty-cycle=0.2,  A/  =2PRF,  1=0,  /=0, 
and  with  Hamming  Window 
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Figure  12.  Plot  of  Cut  along  Frequency  Axis  for  x  =  0  of  the  Cross-ambiguity  Function 
for  the  Stepped  Frequency  Radar ,  Pulse  Number  N=4,  Number  of  Pulses 
Processed  M=4,  Duty  cycle=0.2,  A/=2PRF,  1=0,  1=0,  and  with  Hamming 
Window 


o  ■  t 
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Figure  13.  Plot  of  Cut  along  Time  Axis  for  f=0  of  the  Cross-ambiguity  Function  for  the 
Stepped  Frequency  Radar,  Pulse  Number  N=4,  Number  of  Pulses  Processed 


M-4,  Duty-cycle-0.2  ,A f  -2PRF  ,  1-0,  and  1=0,  and  with  Hamming  Window 
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auto- ambiguity  diagram  for  the  stepped  frequency  radar  without  window  function, 
one  can  find  that  window  functions  which  are  often  used  in  radar  processors  will  decrease 
the  amplitude  of  each  spike  and  spread  the  width  of  it  in  return  for  reducing  the  sidelobes 
of  each  spike.  The  decrease  in  amplitude  will  reduce  the  detection  capability  and  the 
spread  of  the  width  of  the  spike  will  cause  poor  accuracy  and  resolution.  On  the  other 
hand,  reducing  the  sidelobe  of  the  spike  will  enhance  the  detection  of  small  targets  close  to 
the  large  target.  Figure  14  shows  the  bank  of  Doppler  filters  within  a  DFT  processor  . 
This  figure  indicates  that  the  ambiguity  function  for  the  /th  Doppler  filter  is  the  shifted 
version,  along  the  frequency  axis,  of  the  ambiguity  function  for  the  Oth  Doppler  filter.  The 
effect  of  decay  for  different  Doppler  filters  is  caused  by  the  window  function. 

The  improvement  in  clutter  suppression  capability  that  is  achieved  using  pulse 
Doppler  processing,  including  the  use  of  rejecting  I  pulses  for  initialization,  is 
demonstrated  in  Figures  15  thru  28  [Ref.  7].  The  figures  were  drawn  using  the  MATLAB 
computer  program  which  calculates  the  ambiguity  function  defined  by  Equation  47. 
Figures  15  thru  22  show  the  clutter  suppression  capability  for  the  pulse  radar,  and  Figures 
23  thru  28  show  the  clutter  suppression  capability  for  the  stepped  frequency  radar.  In  all 
these  plots,  the  sample  step  of  time  delay  axis  is  4  PRI.  Therefore  only  envelopes  of  the 
ambiguity  diagram  are  shown  in  the  time  delay  axis.  In  the  frequency  shift  axis,  there  are 
many  spikes  occurring  at  the  position  of  multiples  of  PRF,  but  only  a  single  spike  is  shown 
in  the  plots.  The  Hamming  window  is  used  in  all  of  these  plots.  In  Figures  15,  17,  19,  21, 
23,  25,  and  28,  the  number  of  transmitted  pulses  and  processed  pulses  are  equal. 
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frequency  (prf) 

Figure  14.  Plot  of  Cut  along  Frequency  Axis  for  t=0  of  the  Cross-ambiguity  Function 
for  the  Stepped  Frequency  Radar,  Pulse  Number  N=4,  Number  of  Pulses 
Processed  M=4,  Duty  cycle=0.2 ,  A/>2PRF,  1=0,  /=[0, 1,2,3],  and  with 
Hamming  Window 


N=M=100,  and  the  number  of  initial  pulses  rejected,  I,  is  set  to  zero.  In  Figures  16,  18, 
20,  22,  24,  26,  and  28,  the  number  of  transmitted  pulses,  N,  is  increased  to  150,  and  I  is 
set  to  50.  From  Figures  15  thru  18,  in  the  case  of  the  pulse  radar,  the  use  of  initial  pulse 
rejection  leads  to  an  ambiguity  function  which,  along  the  delay  axis,  is  a  cross  correlation 
of  envelope  functions  of  unequal  length.  This  phenomenon  is  not  clear  for  the  stepped 
frequency  radar  due  to  the  different  pulse  frequencies  as  shown  in  Figures  23  and  24. 
Comparing  Figures  19,  21,  25,  and  27  with  Figures  20,  22,  26,  and  28,  respectively,  the 
improvement  in  clutter  suppression  capability  when  utilizing  initial  pulse  rejection  is  clearly 
demonstrated.  Comparison  of  Figures  19  and  25  with  Figures  20  and  26  shows  that  the 
clutter  in  the  0  to  50  PRI  range  will  be  rejected  by  the  use  of  an  initial  pulse  rejection  for  a 
0.1  PRF  Doppler  mismatch.  Figures  21  and  27  show  that,  without  initial  pulse  rejection, 
the  capability  of  Doppler  processing  to  suppress  clutter  is  severely  degraded  when  placing 
30  PRIs  of  delay  (from  the  matched  delay)  when  compared  to  the  cases  shown  in  Figures 
22  and  28.  [Ref.  6] 
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Figure  15.  Cross-ambiguity  Diagram  of  the  Pulse  Radar,  Pulse  Number  N=100, 
Number  of  Pulses  Processed  M=100,  Duty-cycle=0.2,  A/  =0, 1=0, 1=0, 
and  with  Hamming  Window 
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Figure  16.  Cross-ambiguity  Diagram  of  the  Pulse  Radar,  Pulse  Number  N=150, 

Number  of  Pulses  Processed  M=100,  Duty-cycle=0.2,  A/  =0, 1=50, 1=0, 
and  with  Hamming  Window. 


Figure  17.  Plot  of  Cut  along  Time  Axis  for  f=0  of  the  Cross-ambiguity  Function  for  the 
Pulse  Radar,  Pulse  Number  N=100,  Number  of  Pulses  Processed  M=100, 
Duty  cycle=0.2,  A f=  0, 1=0,  /= 0,  and  with  Hamming  Window. 


Pulse  Radar,  Pulse  Number  N=150,  Number  of  Pulses  Processed  M=100, 
Duty-cycle=0.2,  A/  =0, 1=50,  /=0,  and  with  Hamming  Window 
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Figure  19.  Plot  of  Cut  Parallel  to  Time  Axis  for  f=0.02PRF  of  the  Cross-ambiguity 


Function  for  the  Pulse  Radar,  Pulse  Number  N=100,  Number  of  Pulses 
Processed  M=100,  Duty-cycle=0.2,  Af=  0, 1=0,  /= 0,  and  with  Hamming 
Window 


Figure  20.  Plot  of  Cut  Parallel  to  Time  Axis  for  f=0.02PRF  of  the  Cross-ambiguity 


Function  for  the  Pulse  Radar,  Pulse  Number  N=150,  Number  of  Pulses 
Processed  M=100,  Duty-cycle=0.2,  A/=0, 1=50, 1=0,  and  with  Hamming 
Window. 
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Figure  21.  Plot  of  Cut  Parallel  to  Frequency  Axis  for  t  =30PRI  of  the  Cross-ambiguity 
Function  for  thePulse  Radar,  Pulse  Number  N=100,  Number  of  Pulses 
Processed  M=100,  Duty-cycle=0.2,  A/=  0, 1=0, 1=0,  and  with  Hamming 
Window 


Figure  22.  Plot  of  Cut  Parallel  to  Frequency  Axis  for  x  =30PRI  of  the  Cross-ambiguity 
Function  for  the  PulseRadar,  Pulse  Number  N=150,  Number  of  Pulses 
Processed  M=100,  Duty-cycle=0.2,  A f=0, 1=50, 1=0,  and  with  Hamming 
Window. 


Figure  23.  Cross-ambiguity  Diagram  of  the  Stepped  Frequency  Radar,  Pulse  Number 
N=100,  Number  of  Pulses  Processed  M=100,  Duty-cycle=0.2  ,  Af  =0.1PRF, 
1=0,  /=0,  and  with  Hamming  Window 
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Figure  24.  Cross-ambiguity  Diagram  of  the  Stepped  Frequency  Radar,  Pulse  Number 
N=150,  Number  of  Pulses  Processed  M=100,  Duty-cycle=0.2  ,  A/  =0.1PRF, 
1=50, 1=0,  and  with  Hamming  Window 
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Figure  25.  Plot  of  Cut  Parallel  to  Time  Axis  for  f=0.02PRF  of  the  Cross-ambiguity 

Function  for  the  Stepped  Frequency  Radar,  Pulse  Number  N=100,  Number  of 
Pulses  Processed  M=100,  Duty-cycle=0.2 ,  A/’=0.1PRF,  1=0, 1=0,  and  with 
Hamming  Window 


Figure  26.  Plot  of  Cut  Parallel  to  Time  Axis  for  f=0.02PRF  of  the  Cross-ambiguity 


Function  for  the  Stepped  Frequency  Radar,  Pulse  Number  N=150,  Number  of 
Pulses  Processed  M=100,  Duty-cycle=0.2  ,A/’=0.1PRF,  1=50,  /=0,  and  with 
Hamming  Window. 
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Figure  27.  Plot  of  Cut  Parallel  to  Frequency  Axis  for  t=30PRI  of  the  Cross-ambiguity 

Function  for  the  Stepped  Frequency  Radar,  Pulse  Number  N=100,  Number  of 
Pulses  Processed  M=100,  Duty-cycle=0.2  ,  A/'=0.1PRF,  1=0,  7=0,  and  with 
Hamming  Window. 


Figure  28.  Plot  of  Cut  Parallel  to  Frequency  Axis  for  x=30PRI  of  the  Cross-ambiguity 

Function  for  the  Stepped  Frequency  Radar,  Pulse  Number  N=150,  Number  of 
Pulses  Processed  M=100,  Duty-cycle=0.2,  A/  =0.1PRF,  1=50,  7=0,  and  with 
Hamming  Window. 
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IV.  CONCLUSIONS 


This  thesis  investigates  the  ambiguity  functions  of  the  stepped  frequency  waveform 

and  the  stepped  frequency  radar  system.  The  stepped  frequency  waveform  achieves  a  high 

range  resolution  by  coherently  processing  the  returns  from  N  pulses,  each  having  a 

different  earner  frequency  that  changes  by  a  fixed  amount  from  pulse  to  pulse.  The 

mathematical  expression  of  the  auto-ambiguity  function  for  the  stepped  frequency 

waveform  is  derived  from  the  definition.  The  auto-ambiguity  function  is  defined  as  the 

correlation  of  the  complex  envelope  of  the  transmitted  signal  and  the  Doppler  frequency 

shifted  version  of  the  complex  conjugate  of  the  transmitted  signal.  The  3D  plot  of  the 

auto-ambiguity  function  shown  in  Figure  5  is  spiky  like  that  of  a  constant  frequency  pulse 

train  (  shown  in  Figure  A.5)  and  is  tilted  at  an  angle.  Along  the  time  delay  axis,  the 

null-to-null  width  of  each  spike  is  2/N Af  as  compared  to  2/x  for  the  constant  frequency 

pulse  train,  thereby  making  it  possible  to  decrease  the  width  by  increasing  NA/  without 

increasing  its  instantaneous  bandwidth.  However,  for  the  constant  frequency  pulse  train, 

high  resolution  can  be  achieved  by  decreasing  the  pulse  width  and,  thereby,  increasing  the 

instantaneous  bandwidth.  On  the  other  hand,  the  null-to-null  spike  width  of  the  stepped 
frequency  waveform  along  the  frequency  axis  is  equal  to  2/(NxPRI)  which  is  the  same  as 

for  the  constant  frequency  pulse  train.  Thus  frequency  resolution  is  not  improved  by  the 

stepped  frequency  waveform. 
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The  mathematical  expression  of  the  cross-ambiguity  function  for  the  high  PRF 
stepped  frequency  radar  system  including  its  receiver  and  signal  processor  is  also  derived. 
This  mathematical  expression  is  obtained  from  the  output  of  the  Doppler  processor  and  a 
distinct  ambiguity  function  is  defined  for  each  Doppler  processing  filter.  If  weighting  is 
ignored,  the  cross-ambiguity  function  of  the  Oth  Doppler  filter  reduces  to  the 
auto-ambiguity  function  of  the  stepped  frequency  waveform  derived  from  the  traditional 
definition.  Ambiguity  functions  of  other  Doppler  filters  are  shifted  versions  of  the 
auto-ambiguity  function  along  the  frequency  axis.  The  mathematical  expression  has  been 
verified  by  comparing  the  ambiguity  diagram  obtained  by  simulation  which  performs  the 
correlation  method. 

The  improvement  of  the  clutter  suppression  capability  for  the  stepped  frequency 
radar  system  by  rejecting  initial  pulses  is  also  proven.  The  strong  clutter  returns  from  close 
range  can  be  suppressed  by  rejecting  some  initial  pulses,  thereby  reducing  the  interference 
from  the  region.  Four  types  of  basic  waveforms  are  also  discussed  in  the  appendix  for  the 
purpose  of  comparison.  These  are  single  pulse,  periodic  constant  frequency  pulse  train, 
linear  frequency  modulated  pulse,  and  discrete  frequency  modulated  pulse.  The 
mathematical  expressions  of  ambiguity  functions  for  these  waveforms  are  derived  and 
have  been  also  verified  by  simulation. 


46 


APPENDIX  A.  EXAMPLES  OF  THE  AMBIGUITY  FUNCTION  FOR  THE 
BASIC  WAVEFORMS 

In  this  section  the  ambiguity  functions  of  commonly  used  radar  signals  such  as 
single  pulse,  constant  frequency  pulse  train,  linear  frequency  modulated  pulse,  and  discrete 
frequency  modulated  pulse  are  discussed.  The  mathematical  expression  of  the  ambiguity 
function  for  these  waveforms  are  derived  and  verified  by  comparing  with  the  ambiguity 
diagram  obtained  from  the  simulation. 

1.  Single  Pulse  of  the  Sine  Wave 

The  single  pulse  of  the  sine  wave  can  be  defined  as 
St{i)  =  u{t)e^  ,  (A.l) 

where  u(t)  is  the  complex  envelope  of  the  signal,  it  is  defined  as 


u(t)  =  1,  if  0 < t<Ts  , 


=  0 ,  elsewhere,  and 

Ts  is  the  pulse  width. 

(A.2) 

The  ambiguity  function  of  the  single  pulse  can  be  expressed  as 

X(x,fd )  =  J  u(t—  x)u*(t)ej2^d1dt 

— OO 

T 

=|  ej2T?dtdt 

for  0  <  x  <  Ts , 

(A.3-1) 

X(T ,U)  -  C T ,  +  X)  x  X  , 

for  —Ts  <  x  <  0 , 

(A.3-2) 

X(x,fd)  =  0, 

elsewhere. 

(A.3-3) 
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The  absolute  value  of  the  ambiguity  function  of  the  single  pulse  can  be 


expressed  as 


lw*)l=(n-T)x|2gg2 


,  for  0  <  x  <  Ts , 


(A4-1) 


iwai=<r.+t>xi*®ga2i,  for-r,<x<o,  (a«> 

IX(T,^)l=0,  elsewhere.  (A4-3) 

The  ambiguity  function  along  the  Doppler  frequency  shift  axis,  x=0, 

becomes 


\X(0,fd)\=Tsx\^^\ 

(A.  5) 

The  ambiguity  function  along  the  time  delay  axis,/d=0,  becomes 

IX(x,0)l=T,-x  , 

for  0  <x<Ts, 

(A.6-1) 

=  Ts  +  t  , 

for  -Ts  <x<0, 

(A.6-2) 

=  0, 

elsewhere. 

(A.6-3) 

The  ambiguity  diagram,  IX(x/d)l2,  of  the  single  pulse  is  shown  in  Figure 
A.  1 ,  the  contour  plot  is  shown  in  Figure  A.2,  and  plots  of  cut  along  the  frequency  axis  for 
x=0  and  time  delay  axis  for  /d= 0  are  shown  in  Figure  A.3  and  Figure  A.4.  The  plot  for 
zero  time  delay  is  the  square  of  the  spectrum  of  a  rectangular  pulse  which  is  the  square  of 
the  absolute  value  of  the  SINC  function.  The  plot  for  zero  frequency  shift  represents  the 
magnitude  square  of  the  auto-correlation  function  of  a  rectangular  pulse. 


48 


Figure  A.  1.  Ambiguity  Diagram  of  A  Single  Pulse,  Pulse  Width=l. 
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Figure  A.2.  Contour  Plot  of  the  Ambiguity  Diagram  of  A  Single  Pulse,  Pulse  Width=l. 

50 


Magnitude 


frequency  ( 1 / p u 1 s e  -  w i d t h ) 


Figure  A.3 .  Plot  of  Cut  along  Frequency  Axis  for  t=0  of  the  Ambiguity  Function  for 
A  Single  Pulse,  Pulse  Width=l. 


Figure  A.  4.  Plot  of  Cut  along  Time  Axis  for  f=0  of  the  Ambiguity  Function  for 
A  Single  Pulse,  Pulse  Width=l. 
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2.  Constant  Frequency  Pulse  Train 


Constant  frequency  pulse  train  of  N  pulses  can  be  defined  as 


St(t)  =  Z  ]L(t-nTr )eP**t 

n=0 


(A.7) 


where  p(t)  ,  the  complex  envelope  of  the  single  pulse  of  transmitted  signal,  is  defined  as 


where 


\x(t-nTr)  =  1 , 

0, 


if  nTr  ^t<  nTr  +  Ts  , 
elsewhere , 


(A.  8) 


N  :  the  number  of  pulses, 

77  :  pulse  repetition  interval  (PRI), 
fo  :  carrier  frequency. 

The  ambiguity  function  of  periodic  pulse  train  can  be  derived  as 


X(x. ,fd)  =  J  E1 1 ±{t-nTr-%)  E1  ( x*(t-mTr)eJ2^dt 

— o»n=0  m= 0 


(A.9) 


After  changing  variable,  we  obtain  the  ambiguity  function  as 


V-  IN-1  00 

X(zJd)  =  E  E  eW*nTr  J  p(r)p*  n)Tr  -  x)e^“ldt 

n=0m=0  — OO 


(A.  10) 


Now  we  define  the  round  trip  time  in  terms  of  integer  and  fraction  parts 
x  =  (p  +  y)Tr  (A.  11) 


where  p  is  integer  ,y  is  the  fraction,  and  0  <  y  <  1 . 
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After  substituting  Equation  A.  11  in  Equation  A.  10  and  simplifying  ,  we 


obtain  the  expression  for  the  ambiguity  function  as 


X(x,fd)  =  C Ts  -iTr)  x x  eP***'  , 

KfATs^Tr)  n=max  (0>p) 


for  0<y<y , 


(A.12-1) 


(r5  +  (y-  1)7V)  x  ^(A+h-DA)  x  ^„Tr  , 

r)  «=max  (0g^4-l) 


for  1  <y<  1 , 


(A.  12-2) 


elsewhere. 


(A.  12-3) 


We  can  now  obtain  the  ambiguity  function  as  the  following  equation  after 
simplifying  the  summation  term  and  taking  the  absolute  value. 


Iv/-  .r  '*  I  _ /rT'  m  it*  \  . ,  I  wnfdiTs-^Tr)  |  w  |  smvfd(N-\p\)Tr | 

\X(x,fd)\  -  (Ts -iTr)  X  I  X  I  —£rr  I  , 


for  0<y<f, 


(A13-1) 


lw_  Sr  ,  i  \t  \  „ ,  |  smvfd(Ts+(y-l)Tr)  |  w  |  smvfd(N-\p+W)T r  | 

IZ(T,/d)l  -  (Ts  +  (y~  Wr)  X \-^d(TMl_,)Tr)  I  X  I - - 1  , 

for  1  -  y  <  y  <  1 , 


(A  13-2) 


\X(x,fd)\=0, 


elsewhere. 


(A13-3) 
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The  ambiguity  function  along  the  frequency  axis,  i.e.,  T  =  0,  can  be 


expressed  as 


\X(0,fd)\=Tsx\ 


si &WdTs) 
TfdTs 


sin  (nNfdTr) 
sin  (7 jfdTr) 


(A.  14) 


as 


The  ambiguity  function  along  the  time  delay  axis,  i.e.,  fd  =  0,  is  expressed 
\X(x,0)\  =  (N-\p\)(Ts-Yrr), 

for  0  <y<  (A.15-1) 


=  (A/'-l/?  +  ll)((y-l)7V+r*), 

forl-£<y<l,  (A.  15-2) 

=  0  ,  elsewhere.  (A.  15-3) 

The  ambiguity  diagram  of  the  pulse  radar  according  to  Equation  A.  13  is  shown  in 

Figure  A.5  and  the  contour  of  the  ambiguity  diagram  is  shown  in  Figure  A.6.  Plots  of  cut 

along  frequency  axis  for  x=0  and  cut  along  time  delay  axis  for/d=0  are  shown  in  Figure 

A.7  and  A.8.  From  Equation  A.  14  and  Figure  A.7,  spikes  appear  at  frequencies  which  axe 

multiples  of  PRF.  The  first  null  for  each  spike  is  l/NTr  away  from  the  center  of  the  spike. 

Therefore,  each  spike  has  the  null  to  null  width  of  2/AT)  in  the  frequency  axis.  Amplitudes 

of  those  spikes  varies  as  a  SINC  function.  Zero  amplitude  will  appear  at  1/7).  From 

Equation  A.  15  and  Figure  A.8,  it  is  seen  that  there  are  2N-1  spikes  along  the  time  delay 

axis  each  of  which  has  the  width  of  27).  The  plot  of  cut  of  the  ambiguity  function  at  zero 

frequency  shift  represents  the  magnitude  square  of  the  auto-correlation  function  of  a  pulse 

train. 
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Figure  A.5.  Ambiguity  Diagram  of  the  Pulse  Radar,  Number  of  Pulses  N=4, 
Duty-cycle=0.2. 
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8 


Figure  A.6.  Contour  Plot  of  the  Ambiguity  Diagram  for  the  Pulse  Radar,  Number  of 
Pulses  N=4,  Duty-cycle=0.2. 
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Figure  A.7.  Plot  of  Cut  along  Frequency  Axis  for  t=0  of  the  Ambiguity  Function  for  the 
Pulse  Radar,  Number  of  Pulses  N=4,  Duty-cycle=0.2. 


Figure  A.  8.  Plot  of  Cut  along  Time  Axis  for  f=0  of  the  Ambiguity  Function  for  the  Pulse 
Radar,  Number  of  Pulses  N=4,  Duty-cycle=0.2. 
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3.  Linear  Frequency  Modulated  Pulse 

Linear  frequency  modulated  pulse  can  be  mathematically  represented  as 


(A.  16) 


where  u(t)  was  defined  in  Equation  A.  2  and  |i  is  the  rate  of  frequency  change.  The 
ambiguity  function  of  the  linear  frequency  modulated  pulse  can  be  written  as 


X(x,fd)  =  j  u(t- x )ewM2u* 


=  j  uit-xW^e^^dt. 


(A.  17) 


Equation  A.  17  is  exactly  the  same  as  Equation  A.3,  except  that/d  is 


replaced  by/d-|ux.  The  ambiguity  function  becomes 


lw«)l  =  (r,  -x)  X  liggaffi  ■  IotOSkT,, 


(A.18-1) 


for-r^TcO.  (A.  18-2) 


elsewhere. 


(A.  18-3) 
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The  ambiguity  diagram  is  shown  in  Figure  A.9  and  contour  plot  is  shown 
in  Figure  A.  10.  From  Figure  A.2  and  Figure  A.  10,  it  can  be  seen  that  the  ambiguity 
function  of  the  LFM  pulse  is  a  rotated  version  of  the  ambiguity  function  of  a  single  pulse. 
The  angle  of  rotation  depends  upon  the  rate  of  frequency  change  |i  as  shown  in  Figure 
A.11. 

4.  Discrete  Frequency  modulated  Pulse 

The  discrete  frequency  modulated  pulse  consists  of  N  subpulses,  each  at 
different  discrete  carrier  frequency.  Such  a  pulse  can  be  expressed  as 

St(t)  = 2  u(t  -  nT)eJ1^+n^t ,  (A.  1 9) 

n=0 

where  N  is  the  number  of  the  discrete  frequencies  in  the  pulse,  NT  is  the  pulse  width,  and 
u(t)  is  defined  as 

u(t-nT)  =  1,  if  nT<t<(n+l)T, 

=  0,  elsewhere.  (A.20) 

The  envelope  of  the  transmitted  signal  can  be  written  as 

U(t)  =  X  u(t  -  nT)ej2nn£^ .  (A.21) 
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Figure  A.9.  Ambiguity  Diagram  of  the  LFM  Pulse,  Pulse  Width=l. 
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Figure  A.  10.  Contour  Plot  of  the  Ambiguity  Diagram  of  the  LFM  Pulse,  Pulse  Width=l. 
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Figure  A.  1 1  .The  Relation  of  the  Ambiguity  Functions  for  the  Single  Pulse  and  LFM  Pulse. 
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According  to  the  definition  of  the  ambiguity  function,  we  obtain  the 


ambiguity  function  of  the  discrete  frequency  modulated  pulse  as 


X(x ,fd)=j  U{t-z)U*{t)e^ld,  (A.22) 


Now  we  define  the  round  trip  time  in  terms  of  integer  and  fraction  parts, 

x  =  (p+i)T ,  (A.23) 

where  p  is  an  integer  and  y  is  a  fraction  where  0  <  y  <  1 .  We  can  obtain  the  ambiguity 

function  of  the  discrete  frequency  modulated  pulse  after  simplification.  The  expression  of 
the  ambiguity  function  has  the  form  as 


X(x,fd)  =  {(1  -  y)T  x  eMd-p&f)(i+y)TXej2np&fz  x  eMfj-@p+r)Af)W-i+p)T  x 


sinn^-pA/Kl^r  sm7c(/^(2p+7)A/)(A'-l;>l)7'' 
siMt(fd-(2pvf)&f)T  -* 

{y J  x  eMd-ip+i)^lTxj2ti{p+l)&fz  x  gfK(fa-(2p+y+l)^)(N+p)Tx 


^nTUfj-Xp+lWpfT 

K(fd-(p+l)ty)yT 


sin  Ti(fd~(2p+y+l)&f)(N-\ p+l\)T 
sin  it(fcr-(2p+yf  l)A/)r 


}• 


(A.24) 


The  ambiguity  diagram  is  shown  in  Figure  A.  12  and  contour  plot  is  shown 
in  Figure  A.  13.  The  discrete  frequency  modulated  pulse  is  composed  of  the  linear 
frequency  modulated  signal  and  the  pulse  train  with  100%  duty-cycle.  Therefore,  the 
ambiguity  diagram  of  it  has  been  tilted  and  has  many  discrete  spikes. 
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frequency  (PRF) 


Figure  A.  1 2.  Ambiguity  Di 


-1 


time  d 


am  of  the  Discrete  FM  Pulse,  P 
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Figure  A.  13.  Contour  Plot  of  the  Ambiguity  Diagram  of  the  Discrete  FM  Pulse, 
Pulse  Width=l. 


APPENDIX  B.  AMBIGUITY  FUNCTION  PROGRAM  CODES 


%filename :  progl.m 

%  Ambiguity  function  of  the  single  pulse  by  correlation  method 
%update :  27/sep/94 
%by :  Huang,  Jen-chih 

% 

clear; 
tau=l; 
nx=50; 
ny=101; 

t=[0:nx-l]*tau/nx; 
dt=t(2)-t(l); 

z=0; 

fd=linspace(-4/tau,4/tau,ny) 
ul=ones(l,nx); 
for  m=l:ny 

u2=u  1 .  *exp(j  *2*pi*fd(m)  *t); 
c=xcorr(u2,ul).*dt; 
z=[z;(abs(c)).A2]; 
end 

t=[fliplr(-t),t(2:nx)] ; 
figure(l) 
mesh(t,fd,z) 
title(’Ambiguity  diagram  of  single  pulse') 
xlabel('delay  (pulse- width)') 
ylabel('frequency  ( 1/pulse- width)') 
figure(2) 

plot(t,z((ny+l)/2,:))  %plot  of  cut  along  frequency  axis 

xlabel('delay  (pulse- width)') 
ylabel('frequency  ( 1/pulse- width)') 
figure(3) 

plot(fd,z(:,nx))  %plot  of  cut  along  delay  axis 

xlabelfdelay  (pulse- width)') 
ylabel('frequency  ( 1/pulse- width)') 


%  clear  all  the  variable 
%pulse  width  =1 
%#  of  step  in  delay  axis 
%#  of  step  in  frequency  axis 
%time  axis 
%time  step 

%frequency  axis 

%envelope  of  the  transmitted  signal 

%ul  multiply  by  doppler  freq.  shift 
%  correlation  of  ul  &  u2 
%matrix  of  the  ambiguity  function 

%reconstruct  the  time  axis 
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4L  sL»  ml*  «£•  4l#  vL»  «1»  ^  ml*  «1*  «£*  *X>  «!•  «!■  *X»  *X*  «JL»  *sL«  ml*  *sLo  +1*  +1*  ml*  *1*  ^  *1*  «1*  *X»  *1*  vL*  sl^  «£*  *1*  ml*  ml*  ml*  ml*  «X»  *!■  ^  *I»  ml*  ml*  ml*  <Ia  ml*  ml*  ml*  ml*  ml*  ml*  <1#  •!*  >L-  «sL»  *st>  +1*  *L*  •sL*1  4f  4f 

JJm  vfi  yf*  ^  ^  ^  ^  *p  ^X»  <X*  v  v  ^  ^p*  <jw  ^  ^  ^  ^y*  ^  ¥Tw  *T*  ^X*  ^  *y*  *y«  ^  ^  v  v  #TW  v*  ^  *T*  v  *p  •T*  rT*  *7*  •T%  v*  v  'T*  v  *T*  'J' 


^filename :  prog2.m 

%  Ambiguity  function  of  the  linear  frequency  modulation  (LFM)  pulse 
%by  correlation  method 
%  update :  15/aug/94 
%by :  Huang,  Jen-chih 

% 

% 

clear; 
tau=l; 
mu=4; 
nx=50; 
ny=101; 

t=[0:nx-l]*tau/nx; 
dt=t(2)-t(l); 
z=0; 

fd=linspace(-4/tau,4/tau,ny); 
for  m=l:ny 

ul=exp(j*pi*mu.*t.*t); 
u2=conj(ul).*exp(j*2*pi*fd(m)*t); 
c=conv  (u2,fliplr(u  1 )) .  *dt; 
z=[z;(abs(c)).A2]; 
end 
% 

t=[fliplr(-t),t(2:nx)];  %reconstruct  the  time  axis 

figure(l) 

mesh(t,fd,z)  %mesh  plot  of  ambiguity  function 

title('Ambiguity  function  of  the  single  pulse') 

xlabel('delay  (pulse- width)') 

ylabel('frequency  (1/pulse- width)') 

figure(2) 

contour(t,fd,z,40)  %contour  plot  of  the  ambiguity  function 

title('Contour  plot  of  the  signal  pulse') 
xlabel('delay  (pulse- width)') 
ylabel('frequency  ( 1/pulse- width)') 

% 


%clear  all  the  variables 

%pulse  width=l 

%rate  of  frequency  change 

%#  of  step  in  delay  axis 

%#  of  step  in  frequency  axis 

%time  axis 

%time  step 

%frequency  axis 

%envelope  of  the  transmitted  signal 
%ul  multiply  Doppler  frequency  shift 
%correlation  of  ul  &  u2 

%matrix  of  the  mag.  square  of  the  ambiguity  function 
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%  filename :  prog3.m 

%  Ambiguity  function  of  the  stepped  frequency  radar  by  the  correlation  method 
%update :  15/mar/94 
%by :  Huang,  Jen-chih 
% 


clear; 


tau=0.2; 

t0=0.8; 

df=0; 

% 

% 

fs=40; 

N=4; 

T=N*(tau+tO); 

sn=fs*tau; 

tn=fs*tO; 

Tn=sn+tn; 

P=[ones(l,sn),zeros(l,tn)]; 

t=linspace(0,T,N*Tn); 

nx=101; 

fd=linspace(- 1 0, 1 0,nx); 
dt=t(2)-t(l); 

s=0; 


%pulse  width 

%PRI-pulse- width;  PRI=1 
%step  frequency 
%  if  df=0  then  pulse  radar 
%else  stepped  frequency  radar 
%  sample  rate  (Hz) 

%#  of  pulses 

%total  time  period  for  N  pulses 
%#  of  samples  in  pulse  width 
%#  of  samples  between  pulses 
%  of  samples  in  one  PRI 
%  envelope  for  one  PRI 
%time  axis  for  N  pulses 
%#  of  samples  in  frequency  axis 
%frequency  axis 
%time  step 


m=[l,2,3,4];  %step  frequency  sequence 

for  n=l:N 

s=[s,exp(j*2*pi*(m(n)-l)*df*t(Tn*(n-l)+l:Tn*n)).*P];  %transmitted  signal 
end 


z=D; 

for  L=l:mx 
f=fd(L); 

s  1  =conj  (s).*exp(j  *2*pi*f*t); 
s2=s; 

c=conv  (s  1  ,fliplr(s2)).*dt; 

z=[z;(abs(c)).A2]; 

end 

[l,r]=size(z); 
t=[fliplr(-t),t(2:N*Tn)] ; 
figure(l) 
mesh(t,fd,z) 
xlabel('delay  (PRI)') 
ylabel('ffequency  (PRF)') 


%transmitted  signal  multiply  by  doppler  shift 
%the  transmitted  signal 
%correlation  of  si  &  s2 
%matrix  of  the  ambiguity  function 


%resconstruct  time  axis 
%mesh  plot  of  the  ambiguity  function 
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figure(2) 

plot(t,z((nx+l)/2,:))  %plot  of  cut  along  time  axis 

xlabel('delay  (PRI)’) 

grid 

figure(3) 

plot(fd',z(:,(r+l)/2))  %plot  of  cut  along  frequency  axis 

xlabel('frequency  (PRF)') 

grid 

%  ^  ^  ^ . ^ 
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**  *******  *  *  *  ***  **  *  ***  *  *******  *  ******  **  **  **  5fT»fc^C5fC5fc5fc5fC5fc  ***  5jc 2f=  *  **********  **  ***  * 


%filename :  prog4.m 

%Ambiguity  function  of  the  discrete  frequency  pulse 
%update :  18/max/94 
%by :  Huang,  Jen-chih 

% 

% 

clear  %clear  all  the  variables 


tau=l; 

df=l; 

N=4; 

sn=30; 


ny=101; 

t=[0;N*sn-l]*tau/(N*sn); 

dt=t(2)-t(l); 

s=D; 

m=[l,2,3,4]; 
for  n=l:N 


%pulse  width 
%frequency  step 
%#  of  discrete  frequencies 
%  sample  rate  in  delay  axis 
%#  of  sample  in  frequency  axis 
%  time  delay  axis 
%time  step 

%discrete  frequency  sequence 


s=[s,exp(j*2*pi*(m(n)-l)*df*t(sn*(n-l)+l:sn*n))];  %envelope  of  the  transmitted  signal 
end 


z=[]; 

fd=linspace(-4/tau,4/tau,ny);  %ffequency  axis 

for  1=1  :ny 

sl=conj(s).*exp(j*2*pi*fd(l)*t);  %conjugate  of  s  times  Doppler  shift 

s2=s; 


c=conv(sl,fliplr(s2)).*dt;  %correlation  of  si  &  s2 

z=[z;(abs(c)).A2];  %matrix  of  the  ambiguity  function 

end 


t=[fliplr(-t),t(2:N*sn)];  %reconstruct  the  time  delay  axis 

figure(l) 

mesh(t,fd,z)  %mesh  plot  of  the  ambiguity  function 

title(' Ambiguity  function  of  the  discrete  frequency  pulse') 

xlabel('  time  (pulse  width)') 

ylabel('frequency  ( 1/pulse- width) ') 

figure(2) 

contour(t,fd,z,40)  %contour  plot  of  the  ambiguity  function 

title('contour  plot  of  the  ambiguity  function  of  the  discrete  frequency  pulse') 
xlabel('  time  (pulse  width)’) 
ylabel('frequency  ( 1/pulse- width)1) 

% 


*  *  *  *  *  *  *  *  *  *  *  *  **  **  *  *  *  *  *  * **  ***  *  **  *  *  *  *  **  *  ***  *  **  *  *  *  *  *  *  ***  **  **  ****  **  *  5fc  *  *  *  *  *  **  * 


71 


%filename :  progal.m 

%ambiguity  function  of  the  single  pulse  by  equation 
%update :  20/aug/94 
%by :  Huang,  Jen-chih 

% 

% 

clear; 

Ts=l; 

nx=101; 

tau=linspace(-Ts,T  s,nx); 
ny=51; 

fd=linspace(-4/Ts,4/Ts,ny); 

z=D; 

for  m=l  :ny 
f=fd(m); 

x=[]; 

for  n=l:nx 
t=tau(n); 
if  t>=0 

if  f*(Ts-t)==0 
x=[x,(Ts-t)]; 
else 

x=[x,(Ts-t) .  *  (sin(pi*f *(T  s-t)))  ./(pi*f *  (T  s-t))] ; 
end 
else 

if  f*(Ts+t)==0 
x=[x,(Ts+t)]; 
else 

x=[x,  (T  s+t).  *(sin(pi*f*(T  s+t)))./(pi*f*(T  s+t))] ; 
end 
end 
end 

z=[z;x]; 

end 

z=(abs(z)).A2; 
mesh(tau,fd,z) 
xlabel('delay  (pulse-width)') 
ylabel(’fr equency  ( 1/pulse- width) ') 
contour  (tau,fd,  z,40) 

sX-  •\X*  -sL-  sLf  *sL*  vL*  \L»  -sL*  sL*  <s>  •sL*  *sL-  -sL*  vL-  X#  «X*  <L+  sL*  X#  X*  <sl*  si*  X.  «JU  vL«  «J>  «ju 

*’*•  *T%  ^T*  "T*  ■'T*  "T*  v  v  *!*■  ■'T*  *T*  v  #T*  T*  •T*  •T*  T1  *T“  *7*  v  "T*  *1*  v*  'T*  'T*  "T*  *T*  V  ^  ^  /fi  /f*  *71  <7*  *7*  *7»  ^ 


72 


%filename :  proga2.m 

%  Ambiguity  function  of  the  linear  frequency  modulation  pulse  by  equation 

% 

clear; 

Ts=l; 

mu=4; 

nx=101; 

tau=linspace(-T  s,T  s,nx); 
ny=51; 

fd=linspace(-4/T  s  ,4/T  s,ny); 

z=D; 

for  m=l:ny 
x=[]; 

for  n=l:nx 
t=tau(n); 
f=fd(m)-mu*t; 
if  t>=0 

if  f*(Ts-t)==0 
x=[x,(Ts-t)]; 
else 

x=[x,(T  s-t) .  *  (sin(pi*f *  (T  s-t))) ./ (pi*f *  (T  s-t))] ; 
end 
else 

if  f*(Ts+t)=0 
x=[x,(Ts+t)]; 
else 

x=[x,  (T  s+t).  *  (sin(pi*f*  (T  s+t)))./ (pi*f*  (T  s+t))] ; 
end 
end 
end 

z=[z;x]; 

end 

z=(abs(z)).A2; 

mesh(tau,fd,z) 

xlabel('delay  (pulse- width)') 

ylabel('frequency  ( 1/pulse-width) ') 

contour(tau,fd,z,40) 

% 
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Jje  #  #  %  *  #  &  #  2^  *  5k  *  *  5|c  :{:  ;{:  if:  :fc  ;f:  #  5}c  ^c  jf:  :$:  jf:  if:  :£  sf:  #  ;fc  ^c  2{c  jjc  :$c  jfc  ;J:  jf:  ?Jc  jf:  ;Jc  %  ^c  :Jc  5fc  3fc  5}:  jfc  sfc  ;f:  :}:  ;f:  5f:  $z  :fc  sfc  %:  sfc  jJ;  :fc  :f:  sj:  %  ^c  ;j«  :Jc  ifc  %:  ?f: 


%filename :  proga3.m 

%  Ambiguity  function  of  the  pulse  radar  by  equation 
%  update :  20/aug/94 
%by :  Huang,Jen-chih 
clear; 

Ts=0.2; 

Tr=l; 

d=Ts/Tr; 

N=4; 

nx=161; 

tau=linspace(-N*Tr,N*Tr,nx); 

ny=61; 

fd=linspace(- 10, 10,ny); 

p=floor(tau/Tr); 

r=(tau/Tr)-p; 

z=D; 

for  m=l:ny 
x=[]; 
f=fd(m); 
for  n=l:nx 
t=tau(n); 
pp=p(n); 
rr=r(n); 

if  rr>=0  &  rr<d 
if  f==0 

x=[x,(Ts-rr*Tr)*(N-abs(pp))]; 

else 

x  1  =(sin(pi*f*  (T  s-rr*Tr)))./ (pi*f *  (T  s-rr*Tr)) ; 
if  rem(f*Tr,l)==0 
x2=N-abs(pp); 
else 

x2^sin(pi*f*(N-abs(pp))*Tr))V(sin(pi*f*Tr)); 

end 

x=[x,(Ts-rr*Tr).*x  1  *x2] ; 
end 
else 
if  f=0 

x=[x,(Ts+(rr-l)*Tr)*(N-abs(pp+l))]; 

else 

xl=(  sin  (pi*f*  (T  s+(rr- l)*Tr)))./(pi*f*  (T  s+(rr- 1  )*Tr)); 
if  rem(f*Tr,l)==0 
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x2=N-abs(pp+l); 

else 

x2=(sin(pi*f*(N-abs(pp+l))*Tr))./(sin(pi*f*Tr)); 

end 

x=[x,(Ts+(rr-l)*Tr).*xl*x2]; 

end 

else 

x=[x,0]; 

end 

end 

z=[z;x]; 

end 

z=(abs(z)).A2; 

figure(l) 

mesh(tau,fd,z) 

xlabel('delay  (pulse- width)') 

ylabelCfirequency  (1/pulse- width)') 

figure(2) 

contour(tau,fd,z,40) 

% 

************************************************************************ 
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^filename :  proga4.m 

%  Ambiguity  function  of  the  stepped  frequency  radar  by  equation 
%  update :  20/aug/94 
%by :  Huang,  Jen-chih 
clear; 

Ts=0.2; 

Tr=l; 

df=2; 

d=Ts/Tr; 

N=4; 

nx=161; 

tau=linspace(-N*Tr,N*Tr,nx); 

ny=61; 

fd=linspace(-10,10,ny); 

p=floor(tau/Tr); 

r=(tau/Tr)-p; 

z=D; 

for  m=l:ny 
x=[]; 

f=fd(m); 
for  n=l:nx 
t=tau(n); 
pp=p(n); 
rr=r(n); 
fl=f-pp*df; 
f2=f-(2*pp+rr)*df; 
f3=f-(pp+l)*df; 
f4=f-(2*pp+rr+ 1 )  *df ; 
if  rr>=0  &  iT<d 
if  f  1  ==0 
xl=l; 
else 

xl=(sin(pi*fl*(Ts-rr*Tr)))./(pi*fl*(Ts-rr*Tr)); 

end 

if  rem(f2*Tr,l)==0 
x2=N-abs(pp); 
else 

x2=(sin(pi*f2*(N-abs(pp))*Tr))./(sin(pi*f2*Tr)); 

end 

x=[x,(Ts-rr*Tr).*xl*x2]; 
elseifrr<l  &rr>=(l-d) 
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if  £3=0 
xl=l; 
else 

x  1 =(sin(pi*f3*(T  s+(rr- 1  )*Tr)))./(pi*f3*(T  s+(rr- 1  )*Tr)); 
end 

if  rem(f4*Tr,l)==0 
x2=N-abs(pp+l); 
else 

x2=(sin(pi*f4*(N-abs(pp+l))*Tr))./(sin(pi*f4*Tr)); 

end 

x=[x,  (Ts+(rr- 1 )  *Tr) .  *x  1  *x2] ; 

else 

x=[x,0]; 

end 

end 

z=[z;x]; 

end 

z=(abs(z)).A2; 

figure(l) 

mesh(tau,fd,z) 

xlabel('delay  (pulse- width)') 

ylabel('frequency  ( 1/pulse- width)') 

figure(2) 

contour(tau,fd,z,40) 

% 
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^filename :  proga5.m 

%  Ambiguity  function  of  the  discrete  frequency  modulation  pulse  by  equation 
%  update :  20/aug/94 
%by :  Huang,  Jen-chih 
clear; 

N=4; 

T=l/N; 

df=2; 

d=Ts/Tr; 

nx=121; 

tau=linspace(-N*Tr,N*Tr,nx); 

ny=61; 

fd=linspace(-4,4,ny); 

p=floor(tau/Tr); 

r=(tau/Tr)-p; 

z=D; 

form=l:ny 

x=[]; 

f=fd(m); 
forn=l:nx 
t=tau(n); 
pp=p(n); 
rr=r(n); 
fl=f-pp*df; 
f2=f-(2*pp+rr)*df; 
f3=f-(pp+l)*df; 
f4=f-(2*pp+rr+l)*df; 
if  fl*(Ts-rr*Tr)==0 
xl=l; 
else 

xl=(sin(pi*fl*(Ts-rr*Tr)))./(pi*fl*(Ts-rr*Tr)); 

end 

if  rem(f2*Tr,l)==0 
x2=N-abs(pp); 
else 

x2Ksin(pi*f2*(N-abs(pp))*Tr))V(sin(pi*f2*Tr)); 

end 

yl=(Ts-rT*Tr).*exp(j*pi*fl*(Ts+rr*Tr)).*xl.*exp(j*2*pi*pp*df*t) 

.  *exp  (j  *pi  *f2*(N- 1  +pp)  *Tr) .  *x2; 
if  f3*(Ts+(rr-l)*Tr==0 
x3=l; 
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else 

x3=(sin(pi*f3  *  (T  s+(rr- 1 )  *Tr)))./ (pi*f3*(T  s+(rr- 1  )*Tr)) ; 
end 

if  rem(f4*Tr,l)==0 
x4=N-abs(pp+l); 
else 

x4=(sin(pi*f4*(N-abs(pp+l))*Tr))./(sin(pi*f4*Tr)); 

end 

y2=(Ts+(rr-l)*Tr).*exp0*pi*f3*(Ts+(rr-l)*Tr)).*x3.*exp(j5f:2*pi*(pp+l)*dfi:t) 
.  *exp  (j  *pi*f4*  (N+pp)  *Tr) .  *x4; 

x=[x,yl+y2]; 

end 

z=[z;x]; 

end 

z=(abs(z)).A2; 

figure(l) 

mesh(tau,fd,z) 

xlabel('delay  (pulse- width)') 

ylabel('frequency  ( 1/pulse- width) ') 

figure(2) 

contour(tau,fd,z,40) 

% 
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